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Abstract. The perturbative non-linear (NL) effects in the small-x evolution of the gluon densities 
depend crucially on the infrared (IR) regularization. The IR regulator, R c , is determined by the 
scale of the non-perturbative fluctuations of QCD vacuum. From the instanton models and from the 
lattice R c ~ 0.3 fm. For perturbative gluons with the propagation length R c = 0.26 fm the linear 
BFKL gives a good description of the proton structure function F2(x, Q 2 ) in a wide range of x and 
Q 2 . The NL effects turn out to be rather weak and amount to the 10% correction to F2(x,Q 2 ) for 
x < 10- 5 . Much more pronounced NL effects were found in the non-linear model, described in the 
' i '. literature, with a very soft IR regularization corresponding to the IR cutoff at ~ ^qcd- The latter 

<*S j issue is also commented below. 

P. 

Introduction 

The non-perturbative fluctuations of the QCD vacuum [Q]] restrict the phase space for 
the perturbative (real and virtual) gluons of the BFKL cascade [2] thus introducing a new 
scale: the correlation/propagation radius R c of perturbative gluons. From the fits to lattice 
data on field strength correlators R c ~ 0.3 fm [3]. The effects of finite R c are consistently 
incorporated by the color dipole (CD) BFKL equation [0, HI]. The perturbative gluons 
with short propagation length do not walk to large distances, where they supposedly 
fuse together [6]. The fusion probability appears to be controlled by the dimensionless 
parameter R%/ SB, where B stands for the diffraction cone slope fzLS]- 

In this communication we discuss the BFKL [2], BK [6] phenomenology of DIS 
in presence of finite correlation length R c with particular emphasis on a sharpened 
sensitivity of the non-linear effects to the infrared. 



Keywords: small-x evolution, non-linear effects, infrared regularization 
D. PACS: 13.15.+g 13.60.Hb 



Color screening of BFKL gluons 



The finite correlation length implies the Yukawa screened transverse chromo-electric 
field of the relativistic quark, 



i(p)r^g S (p)K l ( P /R c )p/ Pl (1) 



where p is the q — g separation. The kernel Jtf of the CD BFKL equation for the color 
dipole cross section, 



^<x(£,r) = jr«)a(£,r), £=log(l/x), (2) 

is related to the flux of the Weizsacker - Williams gluons \$(p) — $(p + r)\ 2 radi- 
ated by the relativistic gg-dipole r ^,\$\ . The asymptotic freedom dictates that $(p) 
must be calculated with the running QCD charge gs(p) = \/47zas(p) and cts(p) = 
4n/p \og(C 2 /p 2 A 2 QCD ), where C = 1.5. 



DGLAP ordering of dipole sizes and the infrared 



Eq.©, with the BK non-linearity [g] included, greatly simplifies for the DGLAP 
ordering of dipole sizes, r 2 -C p 2 -C R 2 , 



d^,r) = ^a s (r)r 2 j' 2 



<- Rl c dp : 



2(j(^p)- 
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(3) 



Our definition of the profile function in the impact parameter b- space is 



r(§,r,b) 



exp 



IB 



(4) 



where B is the diffraction cone slope and do e i/dt ~ exp [Bt\. Eq.© implies that the 
unitarity limit for o is SnB. 

The diffraction slope for the forward cone in the dipole-nucleon scattering was pre- 
sented in [9] in a very symmetric form 



B(^r)= 1 -(b 2 )= 1 -r 2 + 1 -R 2 N + 2a^. 



(5) 



The dynamical component of B is given by the last term in Eq. © where a^p is the 
Pomeron trajectory slope evaluated first in ffioh (see also |@]). Here we only cite the 
order of magnitude estimate [9J] 
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(6) 



which clearly shows the connection between the dimensionful 0^, and the non- 
perturbative infrared parameter R c . In Eq. © the gluon-probed radius of the proton is 
a phenomenological parameter to be determined from the experiment. The analysis of 
Ref. Ill gives R 2 N « 12 GeV" 2 . 



The function p 2 c(^,p) ~ (Xs(p)G(x,p), where G is the integrated gluon density, is 
fiat in p 2 and the non-linear term in Eq.© is dominated by p ~ R c : 

Thus, the small parameter R^/SB enters the game. To see it in action a partial solution 
to Eq.© is needed. 



Partial solution to GLR-MQ 



The differential form of Eq.© for G is the GLR-MQ equation UJ] 

^<9 7? G(^,T]) = cG(^T])-a(r ] )G 2 (^r ] ), (7) 

where c = 8C F //3 , rj = log [a s {R c ) / a s {p)], 

a(Tj)=a(0)exp[-Tj-A(e T '-l)], (8) 

a(0) = a s (R c )nR 2 c /AfioB and A = 47t/Poa s {R c )- We solve Eq.© making use of the 
exact solution of 

as a boundary condition. To simplify Eq.© we substitute the steeply falling function 
a(r\) witha(T]) = a(0) exp [— yrj]. The partial solution of the simplified problem is found 
readily: 

r( t n) = exp(A^+ 7 T7) 

t»exp(A£) + const' K) 
Here y = c/A with asymptotic value A = 0.4 and 

a M = Rl_nas(R £ ) 
c SB 2N C 

One can see that the gluon fusion mechanism tames the exponential ^ -growth of G(<^ , r\ ) 
but fails to stop the accumulation of large logarithms, r\ = log(l/ (Xs(r)), at very small 
r 2 <7? 2 . 

For vanishing non-linearity, co — » 0, Eq.® matches the exact solution to the CD 
BFKL equation found in Ref . £3 



G ~ exp(A<^ +777 
In the limit ^ — y 00 the gluon density saturates 
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(12) 

and the corresponding dipole cross section reads 

~as(R c 



r 2 



a(r) = SnB-^-2 



OCs(r) 



(13) 



Regime of the additive quark model 



At large r > R c a sort of the additive quark model is recovered: the (anti)quark of the 
dipole r develops its own perturbative gluonic cloud and the pattern of the gluon fusion 
changes dramatically {^\. From Ref.L8Q it follows that the the non-linear correction to 
the dipole cross section is 

-2 ^i v 2 lnlT > MZ,P)a{Z,r) 



80 ~K 2 J C dp 2 Kf(p/R c 

(14) 



%KB 



SnB 

For R c much smaller than the nucleon size <j(R c ) R 2 . Therefore, the magnitude of 
non-linear effects is controlled, like in the case of small dipoles, by the ratio R 2 /SB. 



Small R c - weak non-linearity 

In Ref.dsl] we solved numerically the BFKL and BK equations with purely perturba- 
tive initial conditions and the IR regularization described above. Our finding is that the 
smallness of the ratio R 2 /SnB makes the non-linear effects rather weak even at the low- 
est Bjorken x available at HERA. The linear BFKL with the running coupling and the 
infrared regulator R c = 0.26 fm gives a good description of the proton structure function 
F2(x, Q 2 ) in a wide range of x and Q 2 [8]. For the smallest available x < 10~ 5 the 10% 
NL correction improves the agreement with data, though. 



Soft IR regularization - fully developed non-linearity 

In Ref. lfbUl the non-linear BK-analysis of HERA data on F2(x, Q 2 ) was presented. 
The infrared cutoff for the purely perturbative BK-kernel is deep in the non-perturbative 
region, 

m = 2/Qs - 4G6V- 1 ~ A£ D , (15) 

and the non-linear effects are absolutely important to tame the rapid growth of the linear 
term in the BK equation. The running strong coupling 0Cs(r) at r = rjR in [14] appears 
to be surprisingly small, 0Cs( r iR) — 0.45, thus assuming an applicability of perturbative 
QCD to arbitrarily large distances ~ ^-qcd- However, it is well known that the non- 
perturbative fields form structures with sizes significantly smaller than Aq^ d and local 

field strength much larger than Aq Cd . Instantons are one of them [Q]]. Direct confirmation 
of this picture comes from the lattice [H]. Therefore, the approach developed in ifbUl may 
lead to a very good BK-description of the HERA data but does not agree with the current 
understanding of what is perturbative and non-perturbative in hadronic physics. 



Summary 



It is not surprising that introducing a small propagation length for perturbative gluons 
pushes the nonlinear effects to very small x. More surprising is that within the linear CD 
BFKL approach this small length, R c = 0.26 fm, results in the correct x-dependence of 
Fl(x, Q 2 ) in a wide range of x and Q 2 . The 10% non-linear BK-correction improves the 
agreement with HERA data at smallest available x 

The non-linear BK-description [14] of HERA data which extends perturbative QCD 
to distances ~ ^qcd contradicts to well established non-perturbative phenomena and 
apparently breaks the hierarchy of scales of soft and hard hadronic physics. 
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